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Heavy baryon/meson ratios A c /D° and Aj,/_B° in relativistic heavy ion collisions are studied in 
the quark coalescence model. For heavy baryons, we include production from coalescence of heavy 
quarks with free light quarks as well as with bounded light diquarks that might exist in the strongly 
coupled quark-gluon plasma produced in these collisions. Including the contribution from decays 
of heavy hadron resonances and also that due to fragmentation of heavy quarks that are left in 
the system after coalescence, the resulting \ C /D° and Ab/B° ratios in midrapidity (|y| < 0.5) from 
central Au+Au collisions at y/s NN = 200 GeV are about a factor of five and ten, respectively, larger 
than those given by the thermal model, and about a factor of ten and twelve, respectively, larger 
than corresponding ratios in the PYTHIA model for pp collisions. These ratios are reduced by a 
factor of about 1.6 if there are no diquarks in the quark-gluon plasma. The transverse momentum 
dependence of the heavy baryon/meson ratios is found to be sensitive to the heavy quark mass, 
with the At/B° ratio being much flatter than the A c /D° ratio. The latter peaks at the transverse 
momentum p T ~ 0.8 GeV but the peak shifts to p T ~ 2 GeV in the absence of diquarks. 

PACS numbers: 25.75.-q, 25.75.Cj, 25.75.Dw, 



I. INTRODUCTION 

Recent measurements on nonphotonic electrons from 
decays of midrapidity heavy-flavored mesons in central 
heavy ion collisions at the Relativistic Heavy Ion Collider 
(RHIC) have shown that the nuclear modification factor 
R\a relative to pp collisions is significantly smaller than 
theoretical predictions based on radiative energy loss of 
heavy quarks in produced quark-gluon plasma [3, [H EE 0] • 
To explain the observed small value of R\ A or large 
suppression of heavy meson production in relativistic 
heavy ion collisions, other mechanisms for energy loss 
of heavy quarks in quark-gluon plasma have been pro- 
posed [a,[g,0,HjflLLQ| . Also, it was suggested that the en- 
hancement of charmed baryon production over charmed 
meson production, i.e., the A c /D° ratio, would suppress 
the R e AA of the electrons from heavy mesons [ill, 
This is due to the fact that enhanced charmed baryon 
production reduces the production of charmed mesons 
and thus their contribution to decay electrons. 

In fact, experiments at RHIC have shown that in cen- 
tral heavy ion collisions there is an enhanced produc- 
tion of midrapidity baryons in the intermediate trans- 
verse momentum region not only in light hadrons but 
also in strange hadrons. This enhancement could be 
described by multi-quark dynamics through quark co- 
alescence or recombination [13, [lj|, baryon junc- 
tion loops [ll| [ij], EH , or long-range coherent fields [lj| . 
Therefore, one might expect enhanced heavy baryon pro- 
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duction in the intermediate transverse momentum region 
as well. Since the A c /D° ratio has not been measured in 
heavy ion collisions, it was assumed in Ref. filf to be the 
same as the measured A/K® ratio or in Ref. [lj] to have 
a Gaussian form with a maximum value peaked at the 
transverse momentum p T = 5 GeV. It was then claimed 
that a value of A c /D° ~ l ; which is about a factor of ten 
larger than that in pp collisions at same energy, could 
result in an additional 20-25% suppression of R\ A apart 
from the suppression due to the collisional energy loss of 
heavy quarks in quark-gluon plasma. However, it is not 
yet known what kind of underlying physics can cause 
such a large enhancement of the A c /D° ratio. 

Recently, enhancement of the A c /D° ratio due to the 
existence of [ud] diquarks in strongly coupled quark-gluon 
plasma (QGP) was suggested by Lee et al. [2Cj. This was 
based on the idea of possible existence of quasi bound 
states of quarks and gluons including diquarks in QGP 
when the temperature is between T c and 4T C [21( , where 
T c is the critical temperature for the quark-gluon plasma 
to hadronic matter transition. Based on the quark coa- 
lescence model, it was shown in Ref. [20| that the A c / D° 
ratio could be enhanced by a factor of 4-8 relative to the 
case without diquarks in QGP, depending on the binding 
energy of the light diquark at T c = 175 MeV. The ques- 
tion whether diquarks can exist as quasi bound states in 
cold quark matter was raised earlier in Ref. [22| . where 
the quark-diquark matter was claimed to be energetically 
more favorable than the free quark matter. The idea of 
enhancement of A c baryon production in heavy-ion col- 
lisions was then suggested in Ref. [23[. Based on the 
assumption that the number of diquarks present in a col- 
lision is mainly determined by that in the initial nuclei, 
the author of Ref. [23|] considered the A c /S c ratio and es- 
timated that it could be enhanced by a factor of as large 
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as 80. Since the S c can not be directly measured in near- 
future heavy-ion experiments, the more phenomenologi- 
cally accessible A c /D° ratio is thus addressed in Ref. [201 ] . 

The idea of diquarks 0, [H, [2(| has been widely 
used in hadron physics for describing the mass spectrum, 
electromagnetic properties, and many other properties 
of hadrons (for a review, see, for example, Ref. [2^]). 
although the properties and even the existence of di- 
quarks are still under debate 28] . Among possible signa- 
tures of diquarks in hadrons, the observed large ratio of 
A c /E c ~ 7 in e + e~ collisions was pointed out in Ref. [29j . 
which implies a significantly larger production of scalar 
diquarks than vector diquarks (3(J. (See also Ref. [3l|.) 
In spite of the lack of concrete consensus on the diquark 
nature, it is well-known that the (isoscalar) scalar light 
[ud] diquark of color anti-triplet is the most probable di- 
quark state. Evidence for such scalar diquarks has been 
confirmed by a recent lattice QCD simulation [32]. (See, 
however, Ref. [33j for other lattice QCD calculation on 
diquarks.) 

In diquark models, the A c is usually considered as a 
system consisting of a charm quark and a scalar light 
[ud] diquark. On the other hand, the structure of E c 
and X* is model-dependent: they may have an axial- 
vector light diquark but they also may contain diquarks 
[Qq] made of one heavy quark (Q) and one light quark 
(q) 0, Hll . In Ref. [30] , it is claimed that the ground 
state baryons favor the [Qq] diquark. This is based on 
the observation that the [Qq]q system gives a A c mass 
closer to that from the three-quark calculation, although 
it has a larger mass than that based on the Q[qq] configu- 
ration. Since the color-spin interaction between quarks is 
inversely proportional to the quark masses, the binding 
energies of [Qq] diquarks are smaller than those of scalar 
light [qq] diquarks, and it is not clear whether such [Qq] 
diquarks can exist, in particular, in QGP. (See Ref. [37| 
for a recent discussion on the diquark picture for heavy 
baryons.) 

In this paper, we calculate the Aq/H° ratios (Q = 
b, c and H° stands for D° or B°) in relativistic heavy 
ion collisions based on the quark coalescence model. We 
consider the two cases of Aq being a pure three-quark 
state and of Aq being made of a heavy quark and one 
scalar light diquark. In addition to the direct formation 
of Aq and H°, we also take into account the contribution 
from decays of heavy hadron resonances. Namely, we 
consider the decay of D* mesons for D° production and 
Eq and T,Q baryon decays for Aq production. As we have 
discussed before, we assume that Sq and £q baryons 
are made of three quarks and that neither vector light 
diquarks nor scalar heavy-light diquark exists in QGP. 

This paper is organized as follows. In the next section, 
we briefly discuss heavy hadron production in pp colli- 
sions. A simple thermal model is then used in Sec. Hill to 
show the role of heavy hadron resonances in the Aq /H° 
ratios. The difference between the ratios for charmed 
hadrons and for bottom hadrons is also discussed. In 
Sec. IIVI the coalescence model used in this work is given 



together with both light and heavy quark distribution 
functions. Also described are the calculational meth- 
ods and the contributions from fragmentation of heavy 
quarks that are not used in coalescence. We then show 
in the same section the results for the relative production 
fractions of heavy hadrons and their p T spectra as well as 
the transverse momentum dependence of the Aq /H ra- 
tios. Section fVl contains the conclusions and discussions. 



II. HEAVY HADRON PRODUCTION IN 
PROTON-PROTON COLLISIONS 

Before studying the production of heavy hadrons in 
heavy ion collisions, we briefly discuss their production 
in pp collisions by using the P YTHIA model [381 ] . For pp 
collisions at </s = 200 GeV, the PYTHIA model gives fol- 
lowing relative production fractions of charmed hadrons 
in midrapidity < 0.5), 1 

/(£>+) ~ 0.196, f(D°) ~ 0.607, 

/(£>+) ~ 0.121, /(A c ) ~ 0.076. (1) 

This leads to the particle number ratios 

m @) (2 , 

\ / pp \ / pp 

The above estimated D°/D + ratio can be understood 
largely by considering resonance decays, namely, decays 
of D* mesons into D mesons. First, the D*° meson de- 
cays into D° meson by 100% [39[ since the decay of D*° 
into D + ir~ is prohibited by energy conservation. 2 Sec- 
ond, the D* + meson can decay either into D + or D°. 
The branching ratios are BR(D*+ -> L»°7r+) w 68% and 
BR(D* + — > D + X) w 32%, where sum over any state 
X is understood. If the production rates are equal for 
D and D* mesons aside from the factor due to different 
spin degeneracies, we then would have 

f D°\ 1 + 3 + 3 x 0.68 _ 6.04 
\D+) pp ~ 1 + 3x0.32 ~1.96 ' [) 

where the spin degeneracy of D* is included. We also 
note that the estimated A c /D° ratio is close to that (ss 
0.16) in the SELEX measurement [1|. 

For bottom hadron production in pp collisions at 
y/s NN = 200 GeV, the relative production fractions in 



1 We have confirmed the results reported in Ref. Il'-I . These 
values are also close to the values quoted by Particle Data 
Group (PDG) [H for e+e" annihilation at *Js = 91 GeV: 
f(D+) ~ 0.21, f(D°) ~ 0.54, /(£>+) ~ 0.16, and /(Ac) ~ 0.093. 

2 Note that m(_D*°) = 2007 MeV, m(D*+) = 2010 MeV, m(D°) = 
1865 MeV, m(D+) = 1869 MeV, while m(7r±) = 139.6 MeV and 
m(7r°) = 135 MeV. 
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midrapidity from the PYTHIA model are 

f(B-) ~ 0.101, f(B°) ~ 0.101, /(S s -) ~ 0.030, 
/(£*-) ~ 0.302, ~ 0.301, /(£*-) - °- 089 > 

/(A b ) ~ 0.076, (4) 

leading to 3 

A b /B° » 0.7, (5) 

which is significantly larger than the A c / D° ratio given 
in Eq. ©. 

III. THE THERMAL MODEL 

In relativistic heavy-ion collisions, the A c /D° and 
Ab / B° ratios are expected to be different from those given 
in Eqs. ([2]) and ([5]) for pp collisions. This can be seen us- 
ing a simple thermal model which assumes that in rela- 
tivistic heavy ion collisions charmed and bottom hadrons 
are produced during hadronization of the quark-gluon 
plasma and that they are in both thermal and chemi- 
cal equilibrium at the phase transition temperature Tq- 
The total charm and bottom numbers are, however, de- 
termined by initial hard scattering as thermal production 
of heavy quarks from created quark-gluon plasma is neg- 
ligible in heavy ion collisions at RHIC [4l|. The number 
of heavy hadrons of mass m inside a fireball of volume V 
at Tc is then given by 



N = lQ ^m 2 TK 2 (m/T c ), 



(6) 



where 7q is the heavy-quark fugacity, g is the degeneracy 
of the particle, and K 2 is the modified Bessel function. 

For charmed hadrons, we would then have 
(D° / D + )o = 1, where the subscript means the 
ratio without resonance contributions, if the D* — > D 
decay is not included. Since 



"D°" 



3m 2 D .K 2 {m D ./T c ) 
m 2 D K 2 (m Do /T c ) 



1.47 



(7) 



at Tc — 175 MeV, including decays of D* mesons changes 
the ratio to 



D+ 



1 + (1 + 0.68) x 1.47 
1 + 0.32 x 1.47 



2.36, 



(8) 



if we assume that the same branching ratios for D* de- 
cay to D as in free space. The D° /D + ratio in heavy-ion 



3 We note that the average fractions of bottom mesons and bot- 
tom baryons in pp annihilation at \/s = 1.5 TeV has been esti- 
mated to be f(bu) = f(bd) = 0.399 and /(baryon) = 0.092 
If wc assume that bq mesons include equal number of B and 
B* mesons apart from the spin degeneracy and the number 
of bottom baryons is mostly due to At, then we would have 
Af,/S° ?s 1.1 in pp collisions. 



collisions is thus about 25% smaller than that in pp col- 
lisions. 

For the A c /D° ratio, we have, by using = 
2285 MeV and m D o = 1865 MeV, 



Ac 
D° 



2m 2 A K 2 {m h jT c ) 
ml ( K 2 (m Do /T c ) 



0.24, 



(9) 



where the factor 2 is the spin degeneracy of A c baryon. 
Since the ratio D* / D is 1.47, inclusion of D* resonance 
decays causes a strong reduction of the A c /D° ratio. On 
the other hand, baryon resonances can also contribute to 
the production of A c baryon through their decays and can 
enhance the A c /D° ratio. The major resonance contribu- 
tion to A c comes from E*(2520) of spin-3/2. In thermal 
model, we have 



S^(2520) 

aI 



2 x 3 x 



m%,K 2 [m s »/Tc) 
m\K 2 (m K jT c ) 



(10) 



where the factor 2 comes from the difference in the spin 
degeneracies of A c and S* , and the factor 3 is due to the 
isospin 1 of £*. Since E*(2520) decays almost 100% into 
A c , inclusion of its contribution gives 180% enhancement 
to the number of A c . Another important contribution 
to A c production is from the decay of E c ground state, 
E c (2455) of spin-1/2, via the S c — > A c tt decay. 4 In ther- 
mal model, its abundance relative to that of A r is 



S c (2455) 

"Ar 



3 x 



m\K 2 (m^jTc) 
m\K 2 (m K jTc) 



= 1.26. 



(11) 



Therefore, the contribution from E c (2455) and £*(2520) 
decays increases the A c number by a factor of 3.1 and 
brings the A c to D° ratio close to the naive expectation 
given in Eq. ([9]) based only on directly produced A c and 



D ! 



Ac 
D° 



A C {1 + SJ(2520)/A C + S c (2455)/A c } 



Ac 



D°(l 
1 + 1. 



1.68D*/D) 
+ 1.26 



the 



1 + 1.68 x 1.47 



0.28. 



(12) 



The contribution from higher mass charmed reso- 
nances is negligibly small in the thermal model. For D 
meson production, the next higher state is £>i(2420) of 
J p = 1 + . The branching ratio of the decay of this me- 
son is not well-known except that it decays into D*tt but 
not into Dir. If we assume that it decays into D° via 
an intermediate D* , i.e., Di — > D*ir — > Dirir, the de- 
nominator of Eq. is then modified by the addition 



4 In the strange quark sector, the small mass difference between A 
and E does not allow the strong decay of E into A, and only E° — > 
A7 is allowed. But in heavy quark sector, the mass difference 
between A c and E c is well above the pion mass. 
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of 1.68 x 3 x 0.06 ~ 0.3, where the factor 0.06 is the 
D 1 (2420)/Z?° ratio, giving thus less than 10% enhance- 
ment to D° production. For higher baryon resonances, 
the next excited state is A c (2593). Its contribution to 
the numerator of Eq. (fT2|) is about 0.2, since there is no 
difference in its spin and isospin degeneracy from that 
of A c . Thus it again gives an effect less than 5% to A c 
production. The final value of the A c /D° ratio then es- 
sentially does not change, i.e., the ratio is modified from 
0.28 to 0.27, which is about twice of that in pp collisions. 

The estimated value of the A c /D° ratio from our 
simple thermal model is close to the value obtained in 
the more sophisticated thermal or statistical model of 
Ref. [12], which gives a charmed baryon to meson ratio 
(cqq)/(cq) < 0.25. It is also consistent with the statistical 
model predictions of Ref. [H|, which has A c /D° ~ 0.2 at 
the RHIC energy. 

The Ab/B° ratio in the thermal model can be esti- 
mated in the same way as for the A c /D° ratio. The 
main difference between B° and D° productions is that 
the bottom quark is much heavier than the charm quark. 
As a result, the mass difference between B* and B is 
much smaller than that between D* and D because of 
the manifestation of heavy quark spin symmetry. In fact, 
the mass difference between D* and D is about 140 MeV, 
while that between B* and B is only [3t| 

m(B*) - m(B) « 5325 MeV - 5279 MeV = 46 MeV. 

(13) 

Therefore, unlike D° production, B* mesons cannot de- 
cay into B mesons by strong interaction. 5 This leads 
to a large suppression of B meson production compared 
with the case for D meson production. With m(A(,) = 
5620 MeV, our simple thermal model then gives 



B 



*+ 



B^ 



0.78, ^ 



A b _ 

B° 



0.32. 



(14) 



It indeed shows that (A{,/B°) is somewhat larger than 
(AcAD°) . 

Unlike B production, the production of Ah is affected 
by resonance decays. This is because the mass differ- 
ence between A^ and Ef, approaches to a finite value 
(~ 195 MeV [HI) in the infinite mass limit, although the 
mass difference between E& of spin-1/2 and EJ of spin- 
3/2 becomes smaller as the heavy quark mass increases. 
By using the recent experimental information for the Ef, 
baryons [45| . 

M (E h ) = 5812 MeV, M(E£) = 5833 MeV, (15) 

the Ah / B° ratio becomes 

Ab_ = f Ab_ 
B° ~ \B° 

(16) 



(1 + 3 x 0.35 + 3 x 2 x 0.31) w 1.25, 



5 The decay of B* into B is caused by electromagnetic interactions, 
B* — » B7, and can thus be distinguished in experiments. 



where the factors 0.35 and 0.31 are the Ef,/Af, and E£/Af, 
ratios, respectively. The Ab/B° ratio is seen to be larger 
than the A c / D° ratio by a factor of about 5. Compared 
to its value in pp collisions, the Ab/B° ratio in thermal 
model for heavy ion collisions is close to a factor of two 
larger, similar to the case for the A c /D° ratio. 

There are other hi gher charmed and bottom resonances 
listed by the PDG [39(. These include heavy baryons 
that contain strange quarks such as Sq, Hq, f2g, and 
their resonances. Although most of corresponding bot- 
tom baryons are yet to be discovered, their existence has 
been predicted by the quark model [46j . Using the masses 
listed by the PDG and from the quark model, we find that 
the yield of these baryons is about 60% of the Aq yield 
in the thermal model and is thus appreciable. Including 
these higher resonances in the thermal model changes the 
baryon/meson ratios to 



(17) 



IV. THE COALESCENCE MODEL 



In the thermal model, the relative abundance of par- 
ticles depends only on their masses. Whether diquarks 
exist in QGP thus does not affect the A c /D° and Af,/_B° 
ratios in the thermal model. This is not the case in the 
quark coalescence model as to be shown in the rest of 
this paper. 



A. Formalism 

As in the coalescence model described in Refs. p^ll47|. 
which has been used in Ref. [48j to study the D meson 
transverse momentum spectrum in relativistic heavy ion 
collisions, we assume that quarks, antiquarks, and glu- 
ons in the produced QGP are uniformly distributed in 
a fire cylinder of volume V. Their momentum distribu- 
tions are taken to be thermal in the transverse direction 
but uniform in midrapidity. For the Wigner functions of 
hadrons, they are expressed in terms of Gaussian func- 
tions. Modeling the effect of transverse radial flow of 
the fire cylinder by an effective temperature, the pro- 
duced heavy meson transverse momentum spectrum is 
then given by 



dN 



dp 



9 m- 



M 



V 



dp x dp, 



exp (-fc 2 <r 2 ) S(p 



M 



dN x dN 2 
'dpi dp 2 

-P1-P2), (18) 



where g M is the statistical factor for colored quark and 
antiquark to form a color neutral meson, e.g., g D0 = 1/36 
and g D , = 1/12 for D° and D*°, respectively, and sim- 
ilarly for B° and B*°. The distribution of heavy quarks 
with transverse momentum p\ and light antiquarks with 
transverse momentum p 2 in the fire cylinder frame are 
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denoted by dN%/dpi and dN 2 /dp 2 , respectively. The p M 
is the transverse momentum of produced heavy meson. 
The relative transverse momentum k between the heavy 
quark and light antiquark is defined as 



k 



(m 2 pi - to^) 



(19) 



where m 1 2 are quark masses, and p[ and p' 2 are heavy 
quark and light antiquark transverse momenta, respec- 
tively, defined in the center-of-mass frame of produced 
heavy meson. The width parameter a is related to the 
harmonic oscillator frequency u> by a = I/^/JIlu, where 
/j = m 1 m 2 /(m 1 + m 2 ) is the reduced mass and can be 
related to the size of produced hadrons. Namely, the 
charge rms radius of produced meson is given by 



<r 2 )ch 



3 1 Qim 2 + Q 2 n 
2 [ilo (m 1 + to 2 ) 2 



(20) 



with Q i being the charge of the i-th. (anti)quark. To 
reproduce the root-mean-squared charge radii 0.43 fm of 
D + and 0.62 fm of B + as predicted by the light-front 
quark model of Ref. [4!|, an oscillator frequency uj = 
0.33 GeV is needed. 

We note that contrary to that for mesons consisting 
of only light quarks, the coalescence formula, Eq. (fT8")l . 
does not require heavy and light quarks to have simi- 
lar momenta to form heavy hadrons. Instead, quarks of 
similar velocities have the most probable chance to form 
hadrons. 

Similarly, the coalescence formula for forming Aq, Eq, 
Eq, Eq, and S1q as well as their resonances from three 
quarks are given by 



dN E 
dp B 



9b- 



(2 ^) 6 (a 1 a 2 
V 2 

-k\a\ - 



dp 1 dp 2 dp : 



dNi dN 2 dN 3 
dp 1 dp 2 dp 3 



x exp i 



fcfof) 8{p B -p x -p 2 -p 3 ), 

(21) 



where the index 3 refers to the heavy quark and indices 
1 and 2 to light quarks, and the statistical factor g B has 
values 1/108, 1/36, and 1/18 for Aq, Eq, and Eg, respec- 
tively, and 1/54, 1/18, 1/108, and 1/36 for Eq (E' Q ), E* Q , 
J7q, and Qq, respectively. The relative transverse mo- 
menta defined in the center of mass frame of produced 
baryon are 



fei 
k 2 



1 



m 1 + m 2 
1 



m l + m 2 + m 



(«i 2 Pi - m^) , 

[m 3 (pi + p' 2 ) - (m x + m 2 )p' 3 ] 



The width parameters a i are a i — 1/ ^Jmuj, where 

(m 1 + m 2 )m 3 



Mi 



M 2 



m i 



(22) 



(23) 



and the oscillator frequency uj is related to the baryon 
charge radius by 



(r 2 )ch 



1 



2ui TO, + TO, + TOo 



-Qi 



TOn 



(24) 



The root-mean-squared charge radii of A c and A^ pre- 
dicted by the quark model have similar values of about 
0.39 fm [H3|, and they are reproduced by the oscillator 
frequencies 0.43 GeV and 0.41 GeV, respectively. These 
values are 25-30% larger than the ones for heavy mesons. 
As we shall show in the next subsection, to convert as 
many heavy quarks of small p T into heavy hadrons via 
coalescence as possible 6 requires, however, smaller oscil- 
lator frequencies or larger heavy hadron radii. 

For production of Aq baryons from the coalescence of 
heavy quarks and [ud] diquarks, we consider Aq baryons 
as made of a heavy quark and a light diquark only. Their 
spectrum can thus be obtained from Eq. (|18p by replacing 
g M with g' B = 1/9. Explicitly, it is written as 



dN E 
dp B 



9 b 



dp x dp 2 



dNi dN 2 



V J rz dp x dp 2 
xexp(-k 2 aj q )6{p B - Pl -p 2 ). (25) 

The width parameter a dq or corresponding oscillator fre- 
quency ujdq is determined by fitting the sizes of heavy 
baryons used in the three-quark model. 

The scalar diquarks can also lead to formation of ex- 
cited states of Aq baryons with f* — i and I , which 
can decay into Aq by emitting pions. We note that the 



mass differences between = 



and 



states arc 



small because of heavy quark spin symmetry, e.g., the 
lowest excited states of A c are A(2593) of = \ and 

A c (2625) of j w = | , and thus have a mass difference of 
only 32 MeV. Although the number of these resonances 
is rather small in the thermal model due to their higher 
masses compared to those of ground state Aq baryons, 
it is not the case in the coalescence model if diquarks 
exist in QGP. Since excited Aq baryons have negative 
parities, the scalar diquark inside them is in a relative 
P-wave state with respect to the heavy quark. The coa- 
lescence formula for the production of such hadrons is 



dN B 
dp R 



9b 3^ 



dp x dp 



dNi dN 2 2 2 



dp 1 dp 2 



exp {-k 2 a 2 dq ) 5{p B - p x - p 2 ), 



(26) 



using the Wigner function given in Refs. [51|,[52j]. In the 
above, k is similarly defined as in Eq. (fT9|) and the width 
parameter a dq is taken to have same value as that for the 
ground state Aq. 



6 This is similar to the assumption introduced in the thermal 
model that all charmed quarks are converted into heavy hadrons. 
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B. Quark distribution functions in QGP 

For the quark distribution functions in quark-gluon 
plasma, we use the thermal distribution for light quarks, 
namely, 



dN * \ 9 i V ( irr \ 

= \^3 TO Texp(-m T /T cff ), 



dp 



(27) 



where g q = 6, m T — yP 2 + m\i and p is the transverse 

momentum. Instead of the local temperature and trans- 
verse flow velocity, we use an effective temperature T e g 
(= 200 MeV) and introduce a normalization factor X to 
fix the total number of quarks as in Ref. [HI]. Following 
Ref. we assume that the phase transition tempera- 
ture is T c = 175 MeV and the volume of the quark-gluon 
plasma formed in rapidity \y\ < 0.5 is V = 1, 000 fm 3 for 
central Au+Au collisions at y/s NN = 200 GeV at RHIC. 
This gives the light quark numbers N u — Nd w 245 and 
the strange quark number N s ss 149 if the constituent 
light quark and strange quark masses are taken to be 300 
MeV and 475 MeV, respectively. The resulting normal- 
ization factors X q for light and strange quarks are then 
0.95 and 0.81, respectively. 

Diquarks are also assumed to be in thermal equilib- 
rium and their distribution dN[ ud ^/dp is then the same 
as given in Eq. (|27[) with g^n — 3. The number of di- 
quarks depends on the diquark mass and is in the range 
of 77 - 44 with m [ud] = 455 - 600 MeV If di- 

quarks exist in QGP, the number of free light quarks 
is reduced by twice the number of diquarks so that the 
total number of light quarks is preserved. This would 
cause a reduction of meson yields and an enhancement 
of baryon yields and thus induces larger Aq/H° ratios. 
Results to be presented below are obtained by using the 
diquark mass = 455 MeV, which is in the range of 

diquark mass estimated in Ref. [55]. The sensitivity of 
these results to the diquark mass will, however, be briefly 
discussed. 

For heavy quark transverse momentum distributions in 
midrapidity, we adopt those parameterized in Ref. [|[ for 
central Au+Au collisions at y/s NN — 200 GeV at RHIC, 



dNc 
dp 

dN b 
dp 



19.2 [1 + (p/6) 2 ] 



(1 + p/3.7) 12 [1 + cxp(0.9 - 2p)} ' 



0.0025 



1+, T6 



exp(-p/1.495), (28) 



with the transverse momentum p in unit of GeV. These 
distributions were obtained from heavy quark p T spectra 
from pp collisions at same energy by the number of bi- 
nary collisions (~ 960) in Au+Au collisions. For bottom 
quarks, it was taken from the upper limit of the uncer- 
tainty band of the pQCD predictions in Ref. (5(| for pp 
collisions, as no experimental measurements are available 
at RHIC. The charm quark p T spectrum in pp collisions 



is, on the other hand, determined from fitting simulta- 
neously measure d p T spectrum of charmed mesons from 
d+Au collisions [57j ] and of electrons from heavy meson 
decays in pp collisions. 

Including heavy quark energy loss as estimated in 
Ref. , the resulting heavy quark transverse momentum 
distribution in midrapidity can be parameterized as 



dN, 



EL 



dp 



dN Q 
dp 



l q(p), 



(29) 



with 



L c = 0.8exp(-p/1.2) + 0.6exp(-p/15), 

L b = 0.36 + 0.84exp(-p/10), (30) 

for charm and bottom quarks, respectively, and p again 
in GeV. The heavy quark numbers are then estimated to 
be N c ~ 9.23 and JV& ~ 0.035. The charm quark number 
used here is about a factor three larger than that used in 
Ref. [2(| based on perturbative QCD calculations. This, 
however, does not affect the Aq /H° ratios. For the heavy 

1.5 GeV and m b = 5.0 GeV. 



Results 



quark masses, we use m c 



To calculate the spectra of H° and Aq produced in rel- 
ativistic heavy ion collisions, we consider two scenarios 
for Aq production. In the first model, we assume that 
the Aq baryons have no diquark structure and there are 
no diquarks in QGP either. Therefore, Aq, Sq, Sq, Sq, 
Sq, Sq, J7q, and £Iq baryons are all formed by three- 
quark coalescence. We further include heavy mesons that 
consist of strange quarks such as D s and B s and their 
resonances. Using the same oscillator frequency for both 
heavy mesons and baryons as shown approximately in 
the quark model, we adjust its value to convert all heavy 
quarks of small p T into heavy hadrons via coalescence 
as in the thermal model. This leads to u> = 0.106 GeV 
for charmed hadrons and u> = 0.059 GeV for bottom 
hadrons. The resulting D + and B + charge radii are 
0.74 fm and 1.44 fm, respectively, which are factors of 
1.7 and 2.3 larger than those predicted by the light-front 
quark model of Ref. [4!| . The required smaller oscillator 
frequencies may be partly due to the change in hadron 
radii in medium [58| and partly due to the incomplete 
treatment of production and decay of hadron resonances 
in the present study. 

The second model assumes the diquark structure in 
Aq and the existence of diquarks in QGP as in Ref. [201 ] . 
In this case, Aq baryons are formed by diquark-heavy- 
quark coalescence in addition to three-quark coales- 
cence. 7 Together with the Sq and £q resonance con- 
tributions as well as those from Sq, Sq, Sq, J7q, and 



7 Strictly speaking, the three quarks should coalesce to form the 
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Charmed hadrons 


Bottom hadrons 


Model 


D° D+ D s A c S c /E' c /n c 


B° (B~) B*° (B*~) B° s Bf A b E b /E' b /Q c 


PYTHIA model 


0.607 0.196 0.121 0.076 


0.101 0.3015 0.030 0.089 0.076 


Thermal model 


0.435 0.205 0.179 0.118 0.063 


0.111 0.229 0.049 0.113 0.096 0.062 


Coalescence model (three-quark model) 
(ground state contribution) 
(fragmentation) 


0.348 0.113 0.113 0.288 0.138 
0.051 0.051 0.0264 0.028 0.066 
0.043 0.014 0.009 0.005 


0.053 0.158 0.027 0.080 0.316 0.155 
0.052 0.155 0.026 0.079 0.032 0.073 
0.001 0.003 0.0003 0.001 0.001 


Coalescence model (diquark model) 
(ground state contribution) 
(diquark contribution) 
(fragmentation) 


0.282 0.091 0.123 0.378 0.126 
0.039 0.039 0.028 0.016 0.059 
0.205 

0.048 0.015 0.010 0.005 


0.044 0.131 0.031 0.092 0.385 0.142 
0.040 0.120 0.029 0.088 0.019 0.066 

0.192 

0.004 0.011 0.002 0.004 0.003 



TABLE I: Relative production fractions of midrapidity heavy hadrons produced in central Au+Au collisions at ^/sjvjv = 200 
GeV. The heavy-strange baryon (Hq, Hq, and Qq) yield is suppressed in the PYTHIA model and thus neglected in this work. 



Qq resonances, which are formed only from three-quark 
coalescence as they do not contain the scalar diquark 
structure, excited Aq(| ) and Aq(| ) baryons are also 
included using the P-wave diquark-heavy-quark coales- 
cence. 8 As in the three-quark coalescence model, the os- 
cillator frequencies or width parameters are determined 
by requiring that all low p T heavy quarks are converted 
to hadrons. Besides the common oscillator frequency u> 
for mesons and three-quark baryons, the additional oscil- 
lator frequency for the heavy-quark and diquark system 
is fixed by reproducing the same Aq radius used in the 
three-quark configuration. Then we have w = 0.09 GeV 
for charmed hadrons and ui = 0.049 GeV for bottom 
hadrons, leading to 0.81 fm and 1.58 fm for the D + 
and B + charged radii, respectively. These values lead 
to uj dq = 0.053 GeV and 0.029 GeV, respectively, for 
charmed and bottom baryons consisting of diquarks. In 
obtaining the contribution of resonances to the ground 
state H° and Aq spectra, we neglect for simplicity the 
small momentum shift in their decays. 



i.e. 



In both models, remaining heavy quarks that are not 
converted to hadrons via coalescence, which mostly have 
high p Tl are converted to heavy hadrons by fragmenta- 
tion. This is achieved by assuming that a heavy quark 
of transverse momentum p T fragments into Aq and H° 
with a ratio similar to that given by the PYTHIA model, 



Aq in diquark structure. Here we simply approximate this pro- 
duction process by three-quark coalescence to form non-diquark 
Aq. This thus gives an upper bound for the yield of A c from 
three independent quarks in QGP as pointed out in Ref. [20ll , 
In general, coalescence of three independent quarks can also 
form negative parity baryons. Since this process has a smaller 
contribution than that due to the production of positive parity 
baryons [59ll . it is not considered in present work. 



^) = 0.05exp[-(p T -4.0) 2 /8.0] +0.1, 



with p T in unit of GeV. 



(31) 



1. Relative production fractions of heavy hadrons and 
Aq/H° ratios 

Our results for the relative production fractions of 
heavy hadrons produced in central Au+Au collisions at 
^/s NN = 200 GeV are presented in Table U together with 
those from the PYTHIA model and the thermal model. It 
is seen that the relative production fraction of D° among 
all produced charmed hadrons in the three-quark coales- 
cence model is 0.35, while it is 0.44 and 0.61 for the ther- 
mal model and the PYTHIA model, respectively. This 
fraction reduces to 0.28 in the diquark model. The rela- 
tive production fraction of B° among all produced bot- 
tom hadrons also shows a similar behavior, namely, it is 
0.101, 0.097, 0.052, and 0.043 in the PYTHIA model, the 
thermal model, the three-quark model, and the diquark 
model, respectively. 

For the Aq/H° ratios, we find that without resonance 
contribution A c /D° = 0.55 in the three-quark model and 
A c /D° = 2.2 in the diquark model. These results are dif- 
ferent from those of Ref. [2(| , namely, the A c / D° is larger 
and the enhancement of this ratio due to the presence of 
diquarks is about 4 while it was about 8 in Ref. [20j. 
These differences come from the different values for the 
oscillator frequencies or width parameters as well as the 
fact that in Ref. [2(| the number of free quarks in the 
quark-gluon plasma is not reduced by the presence of di- 
quarks. For bottom hadrons, we obtain Af,/B° = 0.6 in 
the three-quark model, which increases by a factor of 3.5 
to Ab/B° = 2.1 in the diquark model. 
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Including resonances enhances both Aq and H° pro- 
duction. We find that in the three-quark model, about 
70% of D° and about 90% of Aq are produced through 
resonance decays. The importance of resonances is also 
seen in the diquark model. In particular, the diquark 
contribution through negative parity Aq resonances is 
non-negligible. The number of negative parity Aq reso- 
nance of spin-1/2 turns out to be about half of the ground 
state Aq baryon number. Since Aq(^ ) and Aq(| ) 
resonances are almost degenerate in mass, the number of 
negative parity Aq resonances is about two times that of 
the ground state Aq after taking into account the spin 
degeneracy. This should be contrasted with the thermal 
model which predicts suppressed yields of negative parity 
Aq resonances as a result of their larger masses. As a re- 
sult, we find that about 70% of D° and about 75% of Aq 
are produced through resonance decays in the diquark 
model. 

The contribution from heavy quark fragmentation is 
small, particularly to the ratios of particle numbers. This 
is because we have assumed that most hadrons are pro- 
duced by coalescence of low p T heavy quarks. The frag- 
mentation contribution is thus non-negligible only for 
high p T hadrons, which constitute less than 8% of pro- 
duced total heavy hadrons. 

The final Aq / H° ratios after including fragmentation 
and resonance contributions are 

A c /D° = 0.83, A b /B° = 6.00, (32) 

for the three-quark model, and 

A c /D° = 1.34, A b /B a = 8.79, (33) 

for the diquark model. Therefore, the enhancement fac- 
tor for these ratios due to diquarks in QGP is about 1.6 
for both charmed and bottom hadrons. 

Our results thus indicate that the A c /D° (A b /B°) ratio 
in heavy ion collisions predicted by the coalescence model 
is about a factor of 6.4 (8.6) larger than the predictions 
of the PYTHIA model in Section |TT] and about a factor 
of 3.1 (7.0) larger than that from the thermal model in 
Section IIIII The existence of diquarks in QGP further 
enhances this factor to 10.3 (12.6) and 5.0 (10.2), com- 
pared to those from the PYTHIA and thermal models, 
respectively. 

2. Transverse momentum spectra of heavy hadrons 

The transverse momentum spectra of Aq and H° are 
shown in Figs. [T] and [2] for charmed and bottom hadrons, 
respectively, together with the transverse momentum de- 
pendence of the Aq/H° ratio. Solid lines are results 
for the three-quark model and dashed lines are those for 
the diquark model. Compared with the PYTHIA results 
shown by filled circles, the enhancement of Aq yield and 
the suppression of H° yield are evident. The bottom 
hadron spectra are seen to be harder than the charmed 



1 10 2 



io- 



I 1 1 1 1 1 1 1 I 

— three-quark model 

— ■ diquark model 
t PYTHIA 




1 2 3 4 5 
P T (GeV) 



1 2 3 4 5 6 
P T (GeV) 



FIG. 1: Spectra of (a) A c and (b) D°, and (c) the ratio A c /D° 
in midrapidity < 0.5) for central Au+Au collisions at 
,ys NN — 200 GeV. Solid lines are for the three-quark model 
and dashed lines are for the diquark model. Results from the 
PYTHIA model are shown by filled circles. 
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FIG. 2: Same as Fig.Q]for (a) A(, spectrum, (b) B° spectrum, 
and (c) the A /B° ratio. 



hadron spectra due to the flatter bottom quark spectrum 
than the charm quark spectrum. 

Compared with measured p/tt~ and A/ Kg ratios [6(| 



roll |62| j the obtained Aq / H° ratios show very different 
behavior. Although the A c /D° ratio has a similar shape 
as the p/n~ and A/iGj ratios, which could be approx- 
imated by a Gaussian shape as assumed in Ref. [12j, it 
has a much larger width and is thus much flatter than the 
light baryon/meson ratios. Furthermore, the Ab/B° ratio 
is even flatter than the A c /D° ratio. This behavior can 
also be found in the spectra of Aq and H°. Namely, while 
the spectra of Af, and B° have very similar p T dependence 
(Fig. [2]) , those of A c and D° show somewhat different p T 
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dependence (Fig. [J) . All these behaviors originate from 
the large mass of heavy quarks. In the infinite mass limit, 
the momentum of a heavy hadron is completely carried 
by the heavy quark, and the role of light quarks or di- 
quarks is just to give different production probability de- 
pending on the heavy hadron wavefunction. Therefore, 
the Aq/ H° ratio in the infinite heavy mass limit becomes 
a constant and shows no momentum dependence. In the 
case of the A c /D° ratio, finite mass effects are not negli- 
gible but its dependence on the transverse momentum is 
less strong than in the A/ K s ratio. This is clearly seen 
in pp collisions as shown by the PYTHIA model results 
in Figs. [T] and [2] Although the coalescence model shows 
a p T dependence of the Aq/H° ratio, this dependence 
weakens as the heavy quark mass increases. 

We also find that inclusion of diquarks in QGP causes 
a shift of the peak positions in the Aq / H° ratios to lower 
values of p T . In the three-quark model, the A c /D° ratio 
peaks at p T ~ 2 GeV, which shifts to p T ~ 0.8 GeV in the 
diquark model. The change of the peak position is caused 
by both the reduction of the D° meson spectrum and the 
enhancement of the A c spectrum in the diquark model. 
We find, however, that the latter has a larger effect on the 
peak position in the A c / D° ratio. Since the enhancement 
of A c due to light diquarks is more appreciable in the low 
p T region, the peak of the A c / D° ratio in the diquark 
model thus appears at a lower value of p T compared to 
that in the three-quark model. The A^/B ratio also 
shows a shift in the peak position, i.e., fromp T ~ 3.5 GeV 
in the three-quark model to 2 GeV in the diquark model, 
but it is still much flatter than the A c /D° ratio. 

For a more realistic momentum distribution of heavy 
hadrons in relativistic heavy ion collisions after their 
production from QGP, their production and annihila- 
tion due to collisions of and with surrounding particles 
should be investigated by a transport model like the 
AMPT model [H, [13] • However, cross sections for reac- 
tions involving open charmed hadrons are small [65L l6t| , 
and it has been shown in Ref. [5l[ that the effects of 
multiple scattering would be small with such small cross 
sections. Since this work is at the level of an exploratory 
study, such an elaborated investigation is not yet called 
for. 



V. CONCLUSIONS AND DISCUSSIONS 

In this paper, we have studied the Aq/H° (A c /D° and 
Ab/B°) ratios in heavy ion collisions, particularly the en- 
hancement of these ratios due to coalescence of heavy 
quarks with light quarks as well as diquarks that might 
exist in the produced QGP. We have also considered the 
resonance decay effects which are shown to be impor- 
tant for understanding the particle production ratios in 
pp collisions. Our simple estimate based on the thermal 
model as well as p redictions of more sophisticated sta- 
tistical models (42|, EH show that resonance decays also 
plays a crucial role for heavy hadron production in heavy 



ion collisions. We have, therefore, included the contribu- 
tion from resonances in estimating these ratios. 

We have determined the width parameters of hadron 
Wigner functions used in the coalescence model by re- 
quiring that the majority of heavy hadrons at low p T are 
formed from coalescence, which is similar to that assumed 
in the thermal model. The remaining heavy quarks are 
then converted to heavy hadrons via fragmentation as 
in pp collisions. The resulting A c /D° and At,/B° ratios 
in central Au+Au collisions at ^s NN = 200 GeV are 
found to be 1.34 and 8.79, respectively, which are about 
10.3 and 12.6 larger than the predictions of the PYTHIA 
model, and about a factor of 5.0 and 10.2 larger than 
those of the thermal model. The enhancement of the 
Aq/H° ratio due to diquarks is, however, only about 1.6 
for both charmed and bottom Iradrons, which is smaller 
than that predicted in Ref. 20]. The difference mainly 
arises from the different constraints used here and in 
Ref. [2(| for determining the oscillator frequencies or 
width parameters in the hadron Wigner functions and 
from maintaining the same total light quark numbers in 
both the three-quark model and diquark model in the 
present study. 

The enhancement due to diquarks would have been 
larger if we do not take into account the resonance con- 
tribution. This is because of the enhancement of Aq 
baryons from Sq and Sq decays, whose production is 
not affected by the presence of diquarks in QGP. This 
is different from the production of strange A to which £ 
can not decay by strong interactions. Furthermore, in- 
cluding the resonance contribution increases the Af,/B Q 
ratio much more than the A c /D° ratio because B* can- 
not decay into B by strong interactions as a result of 
heavy quark spin symmetry. Finally, this model leads to 
the reduction of D° and B° mesons by factors of about 
2.2 compared to the PYTHIA model. 

Contrary to the PYTHIA model, we have found that 
the yields of heavy baryons with strange quark(s) are 
not small. Their yield was estimated to be about 60% of 
the Aq yield in the thermal model. In the coalescence 
model, this is about 50% and 35% of the Aq yield in the 
three-quark model and diquark model, respectively. 

We have also compared the transverse momentum de- 
pendence of the Aq/H° ratios with that of measured 
A/ Kg ratio and found that the Aq/H° ratios have a 
weaker dependence on the transverse momentum because 
of the massive quarks inside heavy hadrons. In particu- 
lar, this leads to a very weak dependence of the Ab/B° 
ratio on the transverse momentum. 

We have further found that the Aq/D° ratio peaks 
at p T ~ 2 GeV in the three-quark model and at p T ~ 
0.8 GeV in the diquark model. Therefore, the enhance- 
ment of heavy baryon production over heavy meson 
production due to diquarks can mostly be observed at 
low p T region. We have also estimated that the ratio 
A c /£° ~ 23, which is less than that of Ref. ;23i] by a 
factor of ~ 4. 

These results are based on a diquark mass of 445 MeV 
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or a binding energy of 145 MeV. Because of the thermal 
factor, the number of diquarks decreases as the diquark 
mass increases. This would reduce both A c and A& pro- 
duction and increase that of D° and B°, reducing thus 
the Aq/H° ratios. For example, if the diquark mass is 
550 MeV, the enhancement of the A c /D° ratio due to di- 
quarks would reduce from 1.6 to 1.3. The peak position 
in the ratio changes, however, very little. 

In studying resonance production in the coalescence 
model, we have assumed that the effect due to energy 
mismatch between quarks and produced hadron is small. 
Since the coalescence model can be viewed as formation 
of bound states from interacting particles in the system 
with energy mismatch balanced by other particles in the 
system, this would be reasonable if the energy mismatch 
is small. Otherwise, the correction factor can, in princi- 
ple, be determined by evaluating the transition probabil- 
ity in the presence of other particles. Also, the coales- 
cence model can be further improved by imposing energy 
conservation using Brcit-Wigner type spectral functions 
for the produced hadrons [67} . This would naturally lead 
to suppressed production rate due to the energy mis- 
match between quarks and produced hadron. Improving 
the coalescence model, therefore, deserves further stud- 
ies. 

One may also consider contributions to heavy hadron 
production from other diquark states. It is well-known 



that the color-spin interaction for color-sextet diquark 
of spin-1 is, although not strong, attractive. (See, e.g., 
Ref (HI-) This diquark, of course, cannot form a color- 
singlet baryon with one heavy quark, and color-sextet 
diquarks are thus disfavored in phcnomcnological mod- 
els [3l|. However, color-sextet diquarks in QGP can pro- 
duce color-singlet baryons by first forming color-octet or 
color-decuplet baryons and then neutralizing their colors 
by interacting with quarks and gluons in QGP. There- 
fore, if we assume the existence of (weakly-bound) color- 
sextet vector diquarks in QGP, it may bring in addi- 
tional enhancement of Aq baryon production. Since this 
mechanism is very mo del- dependent, the existence and 
properties of such diquark states in QGP should be un- 
derstood before estimating their contributions to heavy 
hadron production. 
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